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Inferring the equation of state of a metastable hard-sphere fluid from the 
equation of state of a hard-sphere mixture at high densities 
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A possible approximate route to obtain the equation of state of the monodispcrsc hard-sphere system in the 
metastable fluid region from the knowledge of the equation of state of a hard-sphere mixture at high densities 
is discussed. The proposal is illustrated by using recent Monte Carlo simulation data for the pressure of a 
binary mixture. It is further shown to exhibit high internal consistency. 



As is well known, simulation studies of the metastable 
fluid branch of the one-component hard-sphere (HS) sys- 
tem, especially near the putative glass transition, are ex- 
tremely difficult. This is due to the natural tendency 
of the system to crystal formation. On the other hand, 
the introduction of a certain degree of size polydispcr- 
sity hinders the formation of local ordered configurations, 
thus providing access to disordered states. The question 
naturally arises as to whether it is possible to infer the 
equation of state of the underlying metastable pure fluid 
from measurements made on the polydisperse one. The 
aim of this work is to provide a constructive affirmative 
answer to that question using recent simulation data of 
the equation of state of a HS binary mixture at high 
densities i 

One of the great assets of using the distribution func- 
tion approach in the theory of liquids is the ability to link 
the structural properties (radial distribution functions) 
and the thermodynamic properties (internal energy, pres- 
sure, chemical potential, . . . ) of a given system. In the 
case of a multicomponent additive HS mixture such a 
link leads to a particularly simple form of the equation 
of state. Let p be the number density of the mixture and 
x i = Pi! Pi where p, is the number density of species i, be 
the mole fraction of species i. Also, let er.y = \(oi +<Jj), 
where the diameter of a sphere of species i is an = cr;, 
be the distance of separation at contact between the cen- 
ters of two interacting particles, one of species i and the 
other one of species j. In terms of the packing fraction 
i] = (-K/Q)p(a 3 ), where (er") = J^iLi x i a ? an d N is the 
number of components, the compressibility factor of the 
mixture Z = p/pksT (p being the pressure, ks the Boltz- 
mann constant, and T the absolute temperature) is given 
by 



N 



1,7 = 1 



22 ^Xjcr- j g lj (a+) 



(1) 



a ' Electronic mail: landres@unex.esl 

http: / / www.uncx.es / eweb / fisteor / andres / 

b J Electronic mail: lsantos@unex.esl 

http: / / www.unex.es / eweb / fisteor / santos / 

c ' Elect ronic mail: |m alopcz@servi dor.unam.mx| 

http: / /miquiztli.cie.unam.mx/xml/tc/ft/mlh/ 



Here gij(cr^j) are the contact values of the radial dis- 
tribution functions <?y(r), r being the distance. Thus, 
one only needs to specify these contact values in order 
to get the equation of state of the mixture. Unfortu- 
nately, no exact analytical results are available for them 
although many valuable approximate expressions have 
been derived or proposed and they may also be deter- 
mined from computer simulations. Among the analytical 
expressions for the contact values, one should especially 
mention the ones which follow from the solution of the 
Percus-Yevick (PY) equation of additive HS mixtures by 
Lebowita^ given by 
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and the results obtained from the Scaled Particle Theory 
(SPT) £2 namely 
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Since neither the PY nor the SPT lead to particularly 
accurate values when compared with simulation results, 
Boublik^ and, independently, Grundkc and Henderson-^ 
and Lee and Levesque^i proposed an interpolation be- 
tween the PY and the SPT contact values, yielding what 
we will refer to as the BGHLL values: 
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These contact values, when substituted into Eq. (JlJ, 
lead to the widely used and rather accurate Boublik 
Mansoori-Carnahan-Starling-Lcland (BMCSL) equa- 
tion of state^ii^ for HS mixtures. Refinements of 
the BGHLL values have been subsequently introduced, 
among others, by Henderson et al.ji2r— , Matyushov and 
Ladanyi,— and Barrio and Solana^ to cope with the so- 
called colloidal limit of binary mixtures. On a differ- 
ent path, but also having to do with the colloidal limit, 
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Viduna and Smit h 24 i 25 have proposed a method to ob- 
tain contact values of the radial distribution functions of 
HS mixtures from a given equation of state. 

In the search for a more general framework able to deal 
with arbitrary spatial dimension d and number of compo- 
nents N, in our previous work 2 ^— we have introduced 
a methodology that provides the contact values of the 
radial distribution functions of d-dimensional multicom- 
ponent hard-core mixtures. In our approach, we have 
taken into account exact results that these contact val- 
ues should satisfy, namely the limit in which one of the 
species, say i, is made of point particles, i.e., <Ji — > 0, the 
limit when all the species have the same size, {<7fc} — > a, 
the limit stemming from a binary mixture in which one 
of the species (say i = 1) is much larger than the other 
one (i.e., <J\/(J2 oo), but occupies a negligible volume 
(i.e., x\{p\j o?) d — > 0), and the limit when the mixture 
is in contact with a hard wall. We have further assumed 
that, at a given packing fraction 77, the dependence of 
9ij( a tj) 011 the parameters {ofc} and {xk} takes place 



only through the scaled quantity 
That is, we have specifically assumed that 
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where, for a given d, the function Q(rj,z) is universal in 
the sense that it is a common function for all the pairs 
regardless of the composition, size distribution, and 
number of components of the mixture. This 'universal- 
ity' assumption has been successfully tested in different 
contexts for HS mixtures in the stable fluid region (see 
Refs. MM for details). 

As can be seen in Eqs. J2])-(|4]), the corresponding con- 
tact values are consistent with the assumption ([6]) , where 
Zij is given by Eq. ([5]) with d = 3 and <?(?7, z%j) has a lin- 
ear (PY) or quadratic (SPT and BGHLL) dependence on 
Zij . In a more general and flexible approach, the contact 
values are given by22r— 
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being the contact value of the pure HS fluid in the SPT 3 
and g p (j]) being the contact value of the pure HS fluid, 
which remains unspecified until one chooses a particu- 
lar equation of state for this latter system. Note that, 
if one chooses g p (r]) = <?p PT (?7), Eq. (J7J) reduces to Eq. 
(j3|). On the other hand, Eq. ([7]) accommodates any de- 
sired function g p (ij), thus extending the contact value of 
the pure fluid to the contact values of a mixture with 
arbitrary composition, number of components, and size 
distribution. Inserting Eq. ([7]) into Eq. ([T]) one obtains 
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where Z p (rf) = 1 + Arig p (r]) is the compressibility factor 
of the pure HS fluid. It can be checked that insertion of 
the PY, SPT, and Carnahan-Starling^ 2 - (CS) expressions 
for Z p {rf) into Eq. ([12]) yields their respective PY, SPT, 
and BMCSL extensions to mixtures. 

Very recently, Odriozola and Bcrthierfi- in the search 
for a thermodynamic signature of the glass transition in 
HSs, have reported Monte Carlo (MC) simulation values 
of the radial distribution functions and the compressibil- 
ity factor for a HS binary mixture up to a packing fraction 
of 77 ~ 0.63. In view of these new results, one can now 
investigate whether (i) our proposal (Eq. (jT2"j0 . used in 
a different way, allows us to infer the equation of state of 
the pure HS fluid in the metastable region, and (ii) the 
related proposal (Eq. (J7J) still holds for high densities. 

Let us begin by noting that Eq. (fT2"j) expresses ^(77) in 
terms of Z p (rf) and known functions of 77. In turn, one 
may invert Eq. (|12p and express Z p (rj) in terms of Z(rj), 
namely 
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Equation (|13|) allows one to obtain indirectly the equation 
of state of the pure HS fluid for high densities beyond the 
freezing transition from the knowledge of the equation of 
state of a given mixture in the same density region. 

Taking as input the MC data of Odriozola and 
Berthier— for the pressure of the binary mixture they 
studied (xi — x 2 = 5, er 2 /<7i = |), we first computed 
the corresponding compressibility factor ^(77) and then 
substituted the values of this latter on the right-hand side 
of Eq. (fl"3"l) to infer the compressibility factor Z p (rf) of the 
pure HS fluid in the metastable fluid region. The results 
appear as full circles in Fig. [TJ where we have also in- 
cluded the accurate molecular dynamics (MD) simulation 
data reported by Kolafa et al.— in the metastable region. 
Figure [1] incorporates as well the curves corresponding to 
the CS compressibility factor and to a free volume com- 
pressibility factor of the form Zf v (rj) = d/(l — rj/r] c ), with 
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FIG. 1. Compressibility factor Z v (j]) of the HS metastable 
fluid as a function of the packing fraction. Filled circles: re- 
sults from Eq. (|13[) : open diamonds: simulation data from Ref. 
l33l : full line: compressibility factor corresponding to the CS 
equation of state; dashed line: compressibility factor arising 
from free volume considerations. The inset shows the same 
data in a different representation. 



FIG. 2. Contact values at high density of the radial distri- 
bution functions of a binary mixture with x\ = X2 = 5 an d 
02/01 = \- Open symbols, with solid lines drawn to guide 
the eye, are the results of Eq. while the filled symbols are 
simulation data from Ref. [l|: <?n(cr^) (red circles); giiip^} 
(black squares); 312 (ujj) (blue triangles). 



d = 3 and where rj c ~ 0.668 is the only fitting param- 
eter. The inset contains the same data and curves, this 
time in the representation 1/Z p (rj) vs 77. It is interesting 
to point out the high degree of consistency in the region 
0.48 < T] < 0.54 between the MD simulation data for the 
true monodisperse HS system 3 ^ and the values inferred 
via Eq. (|13[) from the MC simulation data for the binary 
mixture reported in Ref. QJ. Furthermore, assuming that 
the present results for Z p (rj) are correct for 77 > 0.54, we 
observe that the CS compressibility factor remains accu- 
rate in the metastable fluid region up to 77 ~ 0.58. On 
the other hand, beyond 77 ~ 0.58 it is the free volume fit 
that gives an adequate description. We note that Odrio- 
zola and Bcrthier— fitted the compressibility factor of the 
mixture in the region 0.58 < 77 < 0.63 also with a free 
volume form, except that they obtained ?7 C = 0.669 and 
adjusted the prefactor to be d! = 2.82 instead of d = 3. 
In view of Eq. (JTSJ), we find that the ratio d'/d = 0.94 
can be explained by the factor multiplying Z p , namely 
(<7 2 ) 3 /((T 3 } 2 ~ 0.93, so that both free volume forms are 
fully compatible. 

Once we have examined the possibility to infer the 
compressibility factor for a HS fluid in the metastable 



fluid region from the high density results of the com- 
pressibility factor of a binary mixture with x\ = X2 = 5, 
02/ci = |, we will check on the consistency of our ap- 
proach by testing the validity of Eq. (JT]) . Since it is rea- 
sonable to expect that an accurate value of g p (rj) should 
lead to accurate values of ffy(c^), in order to assess the 
merits of our proposal we take for the contact value of 
the pure HS fluid the one that comes from the Z p (rj) 
that we obtained from the use of Eq. (fT3"|) . The results 
for all pairs 11, 22, and 12 arc shown in Fig. [3J where 
we have also included the simulation data of Odriozola 
and Bcrthier— for these radial distribution functions for 
comparison. As clearly seen in this figure, the agreement 
between our computed values and those of simulation is 
very good. This lends support to the use of Eq. (|7|) also 
for high densities. 

In summary, in this paper we have proposed a new 
route to derive the equation of state of a metastable HS 
fluid. In contrast with our previous work, which rested 
on the notion that the easier system was the pure HS 
fluid while the complicated one was the HS multicompo- 
nent mixture, here we have used the simulation results for 
a HS binary mixture at high density to infer the equa- 
tion of state of a metastable one-component HS fluid. 
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Although not shown, we have also checked that if, in- 
stead of using the simulation data for the pressure of the 
mixture, we use the simulation data for the contact val- 
ues of the radial distribution functions <7ii(<7~l~), 322(0"^"), 
and 312(0" to compute the compressibility factor of the 
mixture Z(rf), the results that we get for Z p (rj) are prac- 
tically indistinguishable from the ones shown in Fig. [T] 

On a different path, we have also shown that our recipe 
to obtain contact values of the radial distribution func- 
tions of multicomponent HS mixtures from the contact 
value of the radial distribution function of the pure HS 
fluid also works for high densities. Given the fact that 
studying mctastablc HS fluids directly in simulations is 
very difficult, the present approach is offered as a possible 
useful alternative. Of course, one cannot reach definite 
conclusions concerning this alternative on the basis of 
the results of the single binary mixture analyzed in this 
paper. However the consistency in very many aspects 
of our approach as well as its relative simplicity are en- 
couraging. We hope that the above serves to motivate 
more simulation studies like the one by Odriozola and 
Berthicri for other mixtures at high density. Once they 
become available, the usefulness of our proposal may be 
further assessed. 
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